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Abstract 

Background: Silencing of the paternal X chromosome (Xp), a phenomenon known as imprinted X-chromosome 
inactivation (l-XCI), characterises, amongst mouse extraembryonic lineages, the primitive endoderm and the 
extraembryonic endoderm (XEN) stem cells derived from it. 

Results: Using a combination of chromatin immunoprecipitation characterisation of histone modifications and 
single-cell expression studies, we show that whilst the Xp in XEN cells, like the inactive X chromosome in other cell 
types, globally accumulates the repressive histone mark H3K27me3, a large number of Xp genes locally lack 
H3K27me3 and escape from l-XCI. In most cases this escape is specific to the XEN cell lineage. Importantly, the 
degree of escape and the genes concerned remain unchanged upon XEN conversion into visceral endoderm, 
suggesting stringent control of l-XCI in XEN derivatives. Surprisingly, chemical inhibition of EZH2, a member of the 
Polycomb repressive complex 2 (PRC2), and subsequent loss of H3K27me3 on the Xp, do not drastically perturb the 
pattern of silencing of Xp genes in XEN cells. 

Conclusions: The observations that we report here suggest that the maintenance of gene expression profiles of 
the inactive Xp in XEN cells involves a tissue-specific mechanism that acts partly independently of PRC2 catalytic 
activity. 
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Background 

During embryonic preimplantation development of fe- 
male mice, starting at the four-cell stage, the paternal 
X chromosome (Xp) undergoes global epigenetic silen- 
cing associated with the establishment of imprinted X 
chromosome inactivation (l-XCI) [1-5]. In the inner cell 
mass (ICM) of the implanting blastocyst, the Xp is then 
reactivated and a random XCI of the Xp or the maternal 
X chromosome (Xm) occurs de novo during the forma- 
tion of the epiblast that will subsequently give rise to the 
adult tissues [1,6]. In contrast, the extraembryonic line- 
ages of the trophectoderm (TE) and the primitive endo- 
derm (PrE) exhibit l-XCI of the Xp, which is maintained 
afterward in the derived lineages of the placenta and the 
yolk sac, respectively [7,8]. 
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Many studies have concentrated on the characterisation 
of random XCI using the ex vivo model of ICM-derived 
female embryonic stem (ES) cells, the differentiation of 
which is accompanied by the onset of X inactivation. In 
these cells, XCI initiates through the c/s-accumulation of 
the Xist noncoding RNA (ncRNA) on the future inactive 
X (Xi), followed by recruitment of Polycomb repressive 
complexes PRC2 and PRC1 and other epigenetic modifi- 
cations, which, together, result in the progressive estab- 
lishment of an inactive state characterised by its extreme 
stability (for review, see [9-11] and references therein). 

In contrast, l-XCI in extraembryonic tissues has been 
less extensively analysed. Studies of developing embryos 
or trophoblast stem (TS) cells derived from the TE [12] 
have revealed that, similarly to the randomly inactivated 
X, the inactive Xp in the TE lineage is associated with Xist 
ncRNA coating, depletion of active epigenetic marks and 
enrichment for the repressive H4K20mel mark, the 
PRC2-dependent H3K27me3 mark and hypermethylation 
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of CpG islands [3,13-18]. Despite these cumulative regula- 
tory locks ensuring the maintenance of Xp silencing, the 
inactive state in the TE seems to be less stable than that of 
adult somatic tissues because transient reactivation of 
some Xp-linked genes occurs spontaneously in both TS 
and TE cells [18]. As a corollary, a large number of X- 
linked genes are expressed from both X chromosomes in 
female TS cells [13]. Intriguingly, the magnitude and 
extent of this escape from I-XCI increase during TE differ- 
entiation into trophoblast giant cells, as revealed by an ac- 
crued frequency of reactivation of different Xp-linked 
transgenes and by reactivation of endogenous Xp loci 
[3,16,19-22]. This relaxed silencing is further exacerbated 
upon depletion of the PRC2 member EED, indicating that 
PRC2, possibly via its H3K27me3 catalytic activity, plays a 
role in the maintenance of I-XCI in the TE lineage [23,24]. 
Collectively, these results suggest that I-XCI is more plas- 
tic than random XCI and indicate the interest in an in- 
depth analysis of the stability of I-XCI in the PrE and its 
derivatives. 

The PrE originates from the ICM and gives rise, after 
differentiation, to the visceral endoderm (VE) and par- 
ietal endoderm (PE) that line the yolk sac, two tissues 
which maintain an inactive Xp [8]. Extraembryonic 
endoderm (XEN) cells have been derived from the PrE 
and show many of the properties of PrE stem cells, in- 
cluding the ability to self-renew indefinitely ex vivo and 
to contribute in a lineage-appropriate manner in vivo, 
although undifferentiated XEN cells contribute predom- 
inantly to the PE and much less efficiently to the VE 
[25-27]. Whilst the inactive Xp in XEN cells is coated 
by the Xist ncRNA, it has been reported not to accumu- 
late the PRC2 complex and associated H3K27me3 [26]. 
EED-mutant embryos, however, show increased and fre- 
quent reactivation of an Xp-linked green fluorescent 
protein (GFP) transgene in cells of both the VE and the 
PE [24]. X-linked GFP reexpression is also observed 
upon loss of Xist coating in the PE, suggesting that both 
Xist ncRNA and PRC2 activity are involved in the main- 
tenance of Xp silencing in differentiated PrE cells [28] . 

In order to compare the characteristics of I-XCI in the 
PrE to the X-inactivation process occurring in other lin- 
eages, we combined two different approaches: (1) profil- 
ing of active and repressive chromatin features along the 
X chromosomes using both chromatin immunoprecipi- 
tation followed by chip hybridisation (ChlP-chip) and 
high-resolution immunofluorescence followed by fluor- 
escent in situ hybridisation (immuno-FISH) studies and 
(2) an analysis of X-linked transcriptional activities at 
the level of single XEN cells by FISH on RNA (RNA- 
FISH) and reverse transcription followed by quantitative 
polymerase chain reaction (RT-qPCR). We observed that 
the Xp in XEN cells, in contrast to findings previously 
reported by other researchers, was globally enriched in 



H3K27me3 compared to the Xm. Intriguingly, we ob- 
served that the topological distribution of H3K27me3 on 
the inactive Xp territory did not coincide strictly with that 
of Xist ncRNA, suggesting that some Xp loci may harbour 
only low levels of H3K27me3. Indeed, similarly to what is 
observed in TS cells, a large number of X-linked genes 
(approximately 15%) lack H3K27me3 on both X chromo- 
somes in XEN cells. These H3K27me3-depleted genes 
tend to be expressed from the inactive Xp in a significant 
number of cells, indicating that they escape from I-XCI. 
The frequency of escape, which appears to be extremely 
variable from gene to gene, was not significantly modified 
after depletion of H3K27 trimethylation using an EZH2- 
specific inhibitor. In addition, I-XCI profiles appear to 
be tightly maintained upon differentiation of XEN cells 
into VE cells. Our observations indicate that whilst wide- 
spread escape from I-XCI emerges as a common charac- 
teristic of extraembryonic stem cells, this escape from 
I-XCI is associated with a different subset of genes in each 
extraembryonic lineage. 

Results 

Inactive X chromosome territory is enriched in H3K27me3 
in XEN cells 

In this study, we analysed one male cell line (GHP7/7) 
and two independently derived female XEN cell lines 
(GHP7/9 and GHP7/3) [26]. All three cell lines harboured 
morphological features typical of XEN cells and expressed 
PrE-specific lineage markers (Additional file 1). The fe- 
male cell lines carry an Xp of 129Sv.Pgkla origin and an 
Xm of 129Sv origin, allowing the identification of allele 
origin through the extensive polymorphisms located in 
the large Pgkla-derived region surrounding the Xist gene. 

In a previous report, the authors indicated that I-XCI 
in XEN cells, unlike that in TS cells, is not accompanied 
by PRC2 recruitment and accumulation of associated 
H3K27me3 on the inactive Xp and suggested that XEN 
cells might express a distinct set of chromatin remodel- 
lers [26]. To address this question we compared, in XEN 
and in TS cells, the mRNA steady-state levels of 40 
known epigenetic regulators using single-cell RT-qPCR 
analysis (Figure 1A and Additional file 2 A and B). On 
the basis of their expression profiles, ten of the forty epi- 
genetic regulators were able to discriminate between the 
two cell populations. These XEN/TS epigenetic discrimi- 
nators included members of the PRC1 and PRC2 com- 
plexes. Unexpectedly, given the previously reported 
findings [26], mRNAs encoding the EZH2 H3K27 meth- 
yltransferase appeared more abundant, globally, in XEN 
cells than in TS cells. 

This observation prompted us to reassess the distribu- 
tion of H3K27me3 in XEN female cells using three- 
dimensional immuno-RNA-FISH for Xist and H3K27me3, 
followed by a deconvolution step (Figure IB). Xist coating 
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Figure 1 H3K27me3 accumulates on the inactive X chromosome in female XEN cells. (A) Heatmap of single-cell, steady-state RNA levels for 
the ten most stringent epigenetic discriminants for trophoblast stem (TS) cells and extraendoderm stem (XEN) cells (P < 5 x 1CT 4 by F-test). The 
complete data set is available in Additional file 2B. R, Regulator preferentially associated with transcription repression; D, Regulator associated 
with a dual function; it can be involved in transcription repression or activation. n = 72 female XEN cells (GHP7/9). (B) Representative image of 
immunofluorescence followed by fluorescence in situ hybridisation on RNA for H3K27me3 (green) and Xist (red) on female XEN cells (GHP7/3 cell 
line). Similar results were obtained with the GHP7/9 cell line (data not shown). Quantification of fluorescence intensities for Xist and H3K27me3 
across the inactive X domain show that the two domains do not strictly overlap. Maximal projections after deconvolution are shown. Scale 
bar = 5 urn. (C) Pie charts showing the percentage of nuclei exhibiting an accumulation of Xist RNA only (red), coaccumulation of Xist RNA and 
H3K27me3 (yellow) and accumulation of H3K27me3 only (green) in female XEN cells (GHP7/3 cell line). (D) Boxplots showing the distribution of 
volumes occupied by Xist RNA and H3K27me3 on the inactive X in XEN cells (GHP7/3 cell line). The two distributions are significantly different. 
P< 0.05 by Kolmogorov-Smirnov test (n > 50). Vertical bars below and above the box-plots show the minimal and maximal values in the cell 
population respectively. 
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of the Xi was associated with a coenrichment in H3K27me3 
in the majority of XEN cells, similarly to what we ob- 
served in TS cells (Figure 1C and Additional file 3). 
Although the volumes of H3K27me3 enrichment and of 
Xist domains on the Xi were not significantly different 
(Figure ID), H3K27me3 and Xist nuclear foci did not 
show strict superimposability in the majority of XEN nu- 
clei (Figure IB), indicating that certain regions of the Xi, 
coated by Xist, may lack H3K27me3. 

These results suggest that XEN and TS cells have dis- 
tinct epigenetic identities. Notably, XEN cells highly ex- 
press the Ezh2 gene. Consistent with this, the Xi in XEN 
cells is globally enriched in H3K27me3. 

A large number of genes on paternal inactive X lack 
H3K27me3 in XEN cells 

In order to confirm these results and identify the genomic re- 
gions of the inactive Xp showing lower levels of H3K27me3, 
we performed ChlP-chip analysis for this histone mark 
on male and female XEN cells using a high-resolution 
microarray covering the X chromosome and a subpart of 
chromosome 17, which served as a ChlP-chip efficiency 
control (see the Methods section for data treatment and 
quality controls). Comparing male and female profiles 
allowed us to detect differences marking specifically the Xi. 

X-linked genes, but not intergenic intervals, appeared 
significantly more enriched in H3K27me3 in female com- 
pared to male XEN cells, indicating that genes constitute 
the preferential targets of H3K27me3 accumulation on 
the Xi (Figure 2A). In addition, the majority of expressed 
genes, as opposed to silent genes, exhibited higher levels 
of H3K27me3 in female compared to male XEN cells. 
This is suggestive of an accumulation of this modification 
on Xi alleles, although we noticed that some expressed 
genes lacked H3K27me3 in both male and female extra- 
embryonic cells, which is indicative of biallelic H3K27me3 
depletion (Figure 2B). We were able to identify these 
H3K27me3-depleted Xi genes (H3K27me3-low) using un- 
biased /c-means clustering based on the statistical com- 
parison of H3K27me3 percentages along gene bodies in 
male and female XEN cells (see Figure 2C, the Methods 
section and Additional files 2C and 4A). Approximately 
15% of expressed genes exhibited such a profile in a fe- 
male XEN cell line (GHP7/9), as opposed to 19.5% that 
showed similar behaviour in TS cells and only about 7% 
in liver cells (Figure 2B and Additional file 4B and C). 
Amongst these H3K27me3-low genes, 19 were common 
to both XEN and TS female cells, and of those, only 2 
were associated with reduced levels of H3K27me3 in the 
female liver (Additional file 4C). 

We then assessed whether this lack of H3K27me3 on 
specific Xi genes in female XEN cells was associated with a 
gain of active histone mark H3K4me2 using a similar 
ChlP-chip approach. Lowered H3K4me2 levels in female 



XEN cells than in male XEN cells were observed at the 
majority of expressed genes, consistent with an absence of 
active histone marks on Xi alleles (Additional file 5A). In 
contrast, H3K27me3-low Xi genes appeared significandy 
enriched in H3K4me2 compared to H3K27me3-high Xi 
genes (P = 0.037 by Kolmogorov-Smirnov test) (Additional 
file 5B and C). We noted, however, that this effect was not 
observed at all H3K27me3-low genes and that some 
H3K27me3-high genes may also show high levels of 
H3K4me2 in female cells, indicating that the anticorrela- 
tion between H3K4me2 and H3K27me3 is not strict. 

On the basis of these ChlP-chip analyses, we conclude 
that both the active and repressive histone marks 
H3K4me2 and H3K27me3 are distributed in a lineage- 
specific manner on Xi genes. Many of the genes lacking 
H3K27me3 on the Xi of XEN cells specifically showed a 
gain of H3K4me2 compared to Xi genes enriched in 
H3K27me3, indicating that some of these H3K27me3-low 
genes adopted the characteristics of active chromatin. 

Genes lacking H3K27me3 in female XEN cells escape from 
imprinted XCI to various degrees 

In order to address whether these active chromatin fea- 
tures are associated with gene expression from the Xi in 
XEN cells, we performed two types of analyses. First, in 
single XEN cells, we measured by allele-specific RT- 
qPCR, the maternal and paternal expression of five 
genes with high H3K27me3 levels (Atp7a, Atrx, Rlim, 
Apool and Sh3bgrl) and five genes with low H3K27me3 
levels (Ftx, Jpx, Pbdcl (also known as 2610029G23Rik), 
Brwd3 and Rps6ka6), as well as Xist transcripts (see 
Additional file 6A for RT-qPCR controls). As expected 
in cells subjected to I-XCI, Xist expression appeared 
to be restricted to the Xp and genes showing high 
H3K27me3 levels were expressed predominantly from 
the active Xm (Figure 3). All five genes lacking global 
H3K27me3, on the other hand, showed a significant 
level of paternal expression in XEN individual cells, indi- 
cating that these genes escaped from I-XCI (Figure 3). 
We noted that the relative levels of expression from both 
the Xp and the Xm were extremely variable, both amongst 
cells and between escaper loci, leading to a highly hetero- 
geneous XEN cell population. Consistent results were 
obtained with mixed cell populations (Additional file 6B). 

Second, we performed RNA-FISH on 13 different X- 
linked regions: 9 regions showing high levels of H3K27me3 
{Fgdl, Atp7a, Pgkl, Elf4-Aifin, Tafl-Ogt, Flna, Cenpi, Tcfe3 
and Rlim) and 4 genes showing low levels of H3K27me3 
{KdmSc, Ftx, Hdac8 and Kdm6a) (Additional file 2D). Im- 
portantly, some of these regions (namely, Tafl-Ogt, Flna, 
Cenpi, Tcfe3, Hdac8 and Kdm6a) displayed different 
H3K27me3 and H3K4me2 levels in XEN than in TS cells 
(Additional file 2C). The comparison of RNA-FISH expres- 
sion profiles between these two extraembryonic cell types 
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Figure 2 Distribution of H3K27me3 on the X chromosome in 
XEN cells. (A) Boxplots of the distribution of H3K27me3 along 
X-linked intergenic regions and X-linked gene bodies of expressed 
X-linked genes (Exp.) or along genes that are not significantly expressed 
(Not exp.) in male and female extraendoderm stem (XEN) cells (GHP7/7 
and GHP7/9 cell lines). Expression data were extracted from the Gene 
Expression Omnibus database [GSE:15519] [29]. n = 1,054 intergenic 
regions and n = 642 X-linked genes. *P < 0.05 by Kolmogorov-Smirnov 
test. (B) Representative examples of H3K27me3 distribution along 
expressed genes enriched in H3K27me3 in female compared to male 
XEN cells ([K27-high]) and along expressed genes showing low levels 
of H3K27me3 in female and in male XEN cells ([K274ow]) (GHP7/7 and 
GHP7/9 cell lines). mm9 University of California Santa Clara (UCSC) 
Genome Browser screenshots are shown. IP, immunoprecipitated DNA; 
IN, input DNA. (C) Scatterplot of H3K27me3 percentages along the 
body of expressed X-linked genes in female (y-axis) relative to male 
(x-axis) XEN cells (GHP7/7 and GHP7/9 cell lines (see Additional file 4A 
for results in the GHP7/3 cell line and Additional file 2C for the 
complete data set). Each dot represents a single gene and its 
respective percentage of H3K27me3 in male and female XEN cells, 
/r-means clustering was applied, which led to the identification of three 
classes of genes. K27-high genes (red dots) are concentrated in the 
upper left quadrant consistently, with them being depleted in 
H3K27me3 in male cells and enriched in H3K27me3 in female cells. 
K27-low genes (blue dots) are significantly depleted in H3K27me3 in 
both male and female cells. In agreement with H3K27me3 marking 
preferentially the silent state, very few expressed genes are enriched 
in H3K27me3 in both male and female cells ([0] genes; indicated by 
maroon dots). These [0] genes may have resulted from the fact that 
gene expression data [29] and our chromatin immunoprecipitation 
followed by chip hybridisation experiment data were obtained from 
different XEN cell lines. The pie chart underneath the scatterplot 
represents the percentage of expressed genes in each H3K27me3 class. 



allowed us to assess the lineage specificity of Xp gene silen- 
cing in extraembryonic lineages (Figure 4). 

In agreement with our chromatin profile analysis, 
most H3K27me3-high Xi genes in each extraembryonic 
cell type were transcribed predominantly from the Xa 
P&t-lacking X) and not from the Xi pfisf-coated X), indi- 
cating that these genes are subject to I-XCI. H3K27me3- 
low genes, on the other hand, appeared to be transcribed 
from the Xi in a significant percentage of nuclei in both 
cell types, indicating that they escaped from I-XCI 
(Figure 4). Two exceptions were observed. Hdac8, which 
we expected to escape from I-XCI in XEN cells based on 
its H3K27me3 profile, showed no significant expression 
from the Xi in these cells (Figure 4B). This may be re- 
lated to its globally highly reduced expression level (low 
steady-state level of mature transcripts [29]), which may 
introduce some stochasticity into transcription initiation. 
Conversely, the Kdm6a gene is transcribed from the Xi 
in a significant percentage of TS nuclei, despite show- 
ing significant levels (about 47%) of H3K27me3 on the 
Xi in this cell line (Figure 4D). However, the proportion 
of nuclei showing this behaviour was significantly 
smaller than that seen in female XEN cells in which 
Kdm6a was poorly enriched (about 10%) in H3K27me3. 
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Figure 3 RT-qPCR analysis of X-linked gene expression in single female XEN cells. (A) Cumulative histograms showing allelic gene 
expression levels of Xist, Atrx, Jpx, Ftx and Pbdd in the 72 individual extraendoderm stem (XEN) cells (GHP7/9 cell line). AU, Arbitrary unit. 
Logarithmic scale is used for Xist. (B) Cumulative histograms of the percentages of female XEN cells showing expression from the inactive 
paternal X (Xp) (dark blue), from the active maternal X (Xm) (red) or from both X chromosomes (light blue) for each indicated gene. X-linked 
genes are grouped according to their H3K27me3 enrichment in female XEN cells. The asterisks mark significant differences found upon analysis 
of the expression profile of a given gene showing low levels of H3K27me3 [K27-low] compared to any of the genes enriched in H3K27me3 
[K27-high] analysed (P < 0.05 by x 2 test). On the left, the diagram depicts the position of analysed genes along the X chromosome. The Pgk1 a 
polymorphic region is shown in grey. 



Kdm6a is known to robustly escape from XCI in somatic 
tissues [30]. 

As previously described for TS cells [13] and other 
adult cell types [31], we observed that maternal allelic 
expression prevailed over that of the paternal allele in 
XEN cells for all the escapers examined, suggesting that, 
for most XCI escapers, the escape mechanism was not 
fully efficient in relieving the inactivation. Indeed, the 
degree of escape, as estimated both by the number of 



nuclei exhibiting transcriptional activity from the Xi 
detectable by RNA-FISH and by the number of cells 
showing significant paternal expression detectable by 
single-cell RT-qPCR, appeared to be extremely variable, 
not only between cell types (for example, Ftx showed 
different levels of escape in TS cells compared to XEN 
cells) but also from gene to gene. This supports the idea 
that expression from the Xi at loci escaping from I-XCI 
is a locally regulated phenomenon. 
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Figure 4 Fluorescence in situ hybridisation on RNA analysis of X-linked gene expression in XEN compared to TS female cells. (A) 

Representative images of fluorescence in situ hybridisation on RNA (RNA-FISH) analysis in extraendoderm stem (XEN) female cells (GHP7/9 cell 
line). The inactive X chromosome is identified through Xist RNA accumulation (green). Primary transcription at the indicated X-linked gene is 
codetected in red. Scale bar = 5 urn. (B) Cumulative histograms of the percentages of nuclei with the depicted expression pattern are shown. 
X-linked genes are grouped according to their H3K27me3 class. [K27-high], Enriched in H3K27me3 in female compared to male XEN cells; 
[K27-low], Low levels of H3K27me3 in female and in male XEN cells. On the left of the histogram, the diagram shows the localisation of RNA-FISH 
probes along the X chromosome, n > 200. (C) The same types of data shown in (A) are given for female trophoblast stem (TS) cells (F3 cell line). 
(D) The same types of data shown in (B) are given for female TS cells. *P < 0.05 by x 2 test for significant differences with expression profiles in 
XEN cells (n > 200). 



Level and extent of escape from imprinted XCI are not 
drastically modified upon inhibition of H3K27me3 in 
XEN cells 

In order to address whether the molecular mechanism 
controlling X-linked gene silencing and/or escape in 
XEN cells involves exclusively the specific distribution of 
H3K27me3 along the Xi, XEN cells were treated with 
GSK126, a potent and specific inhibitor of EZH2, the 
PRC2 catalytic subunit responsible for H3K27me3 ap- 
position [32]. Importantly, the inhibitor had no visible 
impact on XEN cell morphology or growth rate and did 
not significantly change the expression of XEN-specific 
markers, indicating that GSK126 did not have important 
pleiotropic effects on XEN cell homeostasis (Additional 
file 7A and C). Western blot analysis of H3K27me3 
on GSK126-treated compared to untreated XEN cells 
showed that the EZH2 inhibitor induced a progressive 
reduction in the global levels of H3K27me3, leading 
to barely detectable H3K27me3 levels after 5 days of 



treatment (Additional file 7B). Immunofluorescence ana- 
lysis of H3K27me3, coupled with RNA-FISH for Xist, 
confirmed that most XEN cells (65%) treated with 
GSK126 lacked H3K27me3 accumulation on the Xi 
(Xist-coated X) after 5 days treatment (Figure 5A). 

To assess the effect of GSK126-induced depletion of 
H3K27me3 on the expression of genes on the Xi, we 
analysed the expression of Xist, three genes subjected to 
I-XCI (Atp7a, Chid and HmgnS) and three genes escap- 
ing from I-XCI (Ftx, Pbdcl and Jpx) using allelic RT- 
qPCR on pools of XEN cells treated for 5 days with the 
drug. The expression levels of some of the X-linked 
genes, including Xist, were significantly increased after 
GSK126 treatment. Importantly, this increase affected 
exclusively alleles that were already active prior to any 
treatment — that is, the paternal Xist allele and maternal 
alleles of some other X-linked genes — suggesting that 
the inhibitor treatment leads to further derepression of 
active promoters (Figure 5B). In contrast, no significant 
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Figure 5 X-linked gene expression in female XEN cells treated with the EZH2 inhibitor GSK126. (A) Representative images produced by 
immunostaining followed by fluorescence in situ hybridisation on RNA for H3K27me3 (green) and Xist (red) on female extraendoderm stem (XEN) 
cells (GHP7/9 cell line) treated or not with 2 uM GSK126 for 5 days. Arrowheads indicate the nuclear position of the inactive X chromosome 
coated with Xist RNA. The percentage of visible H3K27me3 accumulation on the X/sf-coated X chromosome in each condition is indicated, n > 50. 
Scale bar = 5 um. (B) Cumulative histograms showing allelic gene expression levels of Xist, Atp7a, Chici, Hmgn5, Jpx, Ftx and Pbdd in female XEN 
cells (GHP7/9 cell line) treated (+) or not (-) with 2 uM GSK1 26 for 5 days. Mean values and standard deviations of two independent experiments 
are shown for each gene. Values have been standardised by the ubiquitously expressed RplpO gene. On the left, the diagram shows the position 
of the genes analysed by reverse transcription followed by quantitative polymerase chain reaction along the X chromosome. No significant 
difference in the relative levels of paternal X (Xp) chromosome expression could be detected between treated and untreated cells (P > 0.05 by x 2 
test). AU, Arbitrary unit; [K27-high], Enriched in H3K27me3 in female compared to male XEN cells; [K27-low], Low levels of H3K27me3 in female 
and in male XEN cells; Xm, Maternal X chromosome. 



change in expression levels of Xp alleles (except from 
Xist, as mentioned previously) was observed in GSK126- 
treated compared to untreated XEN cells, indicating that 
the marked reduction of H3K27me3 on the Xi of XEN 
cells did not significantly modify the efficiency of Xi 
gene silencing (Figure 5B). Similar results were obtained 
with two other H3K27me3-low genes (Kdm6a and 
Nkap) at the single-cell level using RNA-FISH following 
H3K27me3 immunostaining (Additional file 7D and E). 
Note that neither of the two analyses completely ex- 
cludes the possibility that highly localised H3K27me3 
may have remained associated with gene promoters, 
after GSK126 treatment, in some cells. 



This analysis shows clearly that, although H3K27me3 
enrichment allows the selective identification/segregation 
of genes subjected to I-XCI compared to genes escaping 
from I-XCI, this histone mark is probably not, by itself, 
determinant for X-linked gene expression control. This 
implies that other epigenetic marks or transcriptional 
mechanisms might be implicated in the regulation of Xi 
expression in XEN cells. 

Escape from imprinted X-inactivation is not modified 
upon XEN cell conversion into visceral endoderm 

Because observations of expression of genes from the 
Xp have previously been reported in the VE of EED 
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mutants [24,33], we hypothesized that this extraembry- 
onic tissue might be prone to labile Xp expression. Thus, 
we treated our female XEN cells with bone morpho- 
genetic protein 4 (BMP4), a factor known to trigger 
XEN differentiation towards the VE fate [34,35], and 
then assessed the status of transcription at the Xi. After 
5 days of induction, XEN cells formed a homogeneous, 
epithelium-like layer of cells characteristic of VE conver- 
sion (Figure 6A). Our expression analysis of lineage de- 
terminants at the single-cell level shows that the vast 



majority of cells had repressed PrE- and PE-specific 
markers and activated VE-specific genes, thus confirming 
the VE differentiation (Figures 6B and C and Additional 
file 2E). Transcription from the Xi, which we assessed 
by RNA-FISH, was not significantly induced upon VE 
conversion at the H3K27me3-high genes Tafl, Rlim and 
Pgkl or at the H3K27me3-low gene KdmSc (Figures 6D 
and E). This observation indicates that, unlike what we 
saw during TS cell differentiation (Additional files 2E 
and 8), no change in Xi expression — at least for the 
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Figure 6 X-linked gene expression upon conversion of XEN cells into visceral endoderm. (A) Representative photographs illustrating the 
epithelium-like cell morphology observed after 5 days of bone morphogenetic protein 4 (BMP4) induction of extraendoderm stem (XEN) female 
cells (GHP7/9 cell line) are shown. (B) Heatmap of single-cell, steady-state RNA levels for the 1 8 most stringent discriminants of cells converted 
into visceral endoderm and undifferentiated XEN cells is shown (P < 1CT 3 by F-test). Hierarchical clustering has been applied. The complete data 
set is available in Additional file 2E. n = 72 female XEN cells (red) and n = 82 XEN cells treated with BMP4 (25 uM) for 5 days (yellow). The expected 
tissue specificity of each gene is indicated: EPI, Epiblast; ICM, Inner cell mass; PE, Parietal endoderm; PI. Lab, Placental labyrinth; PrE, Primitive 
endoderm; TE, Trophectoderm; TGC, Trophoblast giant cell; VE, Visceral endoderm. The colour scale is the same as the one shown in Figure 1 A. 
(C) Three-dimensional projections of principal components (PCs) 1, 2 and 3 of single-cell expression profiles of genes shown in (B). (D) Two-colour 
fluorescence in situ hybridisation on RNA analysis for the indicated genes in XEN cells treated with BMP4 for 5 days. Arrowheads point to transcription 
signals associated with the X/sf-coated inactive X chromosome. Scale bar = 5 um. (E) Cumulative histogram showing the percentage of nuclei with 
the depicted expression pattern in XEN cells treated with BMP4 for 5 days (D5). Results obtained in untreated XEN cells (DO) are also shown to facilitate 
the comparison. No significant differences between X-linked gene expression profiles in BMP4 treated compared to untreated XEN cells could be 
detected (x 2 test), n > 100. [K27-high], Enriched in H3K27me3 in female compared to male XEN cells; [K27-low], Low levels of H3K27me3 in female 
and in male XEN cells; Not exp., genes not significantly expressed; Xa, Active X; Xi, Inactive X. 
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genes we tested — is associated with XEN differentiation 
into VE. 

Discussion 

Until now, it has been unclear whether differences in the 
silencing of genes on the Xi in the extraembryonic TE 
lineage, compared to tissues such as liver or embryonic 
kidney [31], are linked to the implementation of differ- 
ent mechanisms of X inactivation (that is, imprinted ver- 
sus random) or whether they are associated with tissue 
specificity. In this study, we identified a large number 
of genes that escape from I-XCI solely in female XEN 
cells and not in female TS cells. Unlike the relaxation 
of I-XCI observed upon TS/TE cell differentiation 
(Additional file 8 and [16]), we observed that gene ex- 
pression profiles from the inactive Xp in XEN cells 
appeared to be tightly maintained upon XEN conversion 
into VE cells in vitro. Hadjantonakis et al. made the 
same observation for yolk sac formation in vivo [19]. 
These observations strongly support the existence of dis- 
tinct regulation of Xp silencing in each extraembryonic 
lineage of the blastocyst. According to our single-cell 
RT-qPCR analysis, the expression of I-XCI escapers does 
not correlate or anticorrelate significantly with the 
expression of any of the lineage biomarkers we tested. 
This is true in both XEN and TS cells, suggesting either 
that other, as yet unidentified regulators must be in- 
volved and/or that tissue specificity is established after 
higher-level regulation of these factors (Additional file 
2B). For instance, such tissue specificity could be medi- 
ated by differences in the topological distribution of 
some of these factors along the inactive Xp in the PrE 
compared to the TE. Specific clustering of I-XCI es- 
capers along the X chromosome was not observed in 
XEN cells, which is suggestive of gene-centred regula- 
tion (Additional file 2C). Gene-specific, ds-acting re- 
gions, either promoting escape from XCI or preventing 
the spreading of active chromatin into adjacent silent 
domain, have been identified in cells subject to random 
XCI, suggesting that similar ds-acting sequences may be 
recruited in extraembryonic tissues [36-38]. Adjacent 
ncRNAs and/or CTCF binding (CCCTC-binding factor 
(zinc finger protein)) may participate in insulating XCI 
escapers from globally inactive chromatin [39-44]. 

Another important observation concerns the variation 
in the degree of escape from gene to gene and, for cer- 
tain genes, between XEN cell lines. Such variability has 
been reported, albeit to a lesser extent, in other tissues 
and organisms [31,45-48]. This characteristic may be 
symptomatic of an accrued lability in the silencing of 
discrete genomic locations and/or a manifestation of 
dynamic transitions between on/off states at paternal 
alleles, as observed previously in TS cells [18]. Ontology 
analyses of XEN-specific escapers did not reveal any 



striking overrepresentation of particular classes of bio- 
logical function compared to that of genes stably sub- 
jected to I-XCI. Neither were we able to identify genes 
known to be involved specifically in yolk sac metabolism 
or development amongst the XEN-specific escapers. 
This could support the hypothesis that genes that are 
not crucial to XEN cell homeostasis are able to tolerate 
low levels of biallelic expression, whereas genes critical 
for XEN function might require tight control of their 
expression level. More detailed gene-by-gene functional 
studies are needed to definitively address this point 
and the role of such large-scale escape in PrE lineage 
development. 

Our results show that the repressive histone mark 
H3K27me3 accumulates on the Xi in XEN cells in a 
manner similar to that observed in other cellular con- 
texts, notably in its preferential targeting to the body 
of inactivated genes and in its exclusion from XCI es- 
capers [17,29,31,49,50]. Strikingly, catalytic inactivation 
of EZH2 does not significantly affect the silencing of Xp 
genes in XEN cells, suggesting that H3K27me3 is not— 
at least not on its own — responsible for the maintenance 
of Xp gene silencing in XEN cells. Rather, the low levels 
of H3K27me3 at I-XCI escapers may be a by-product of 
the initial definition of loci escaping from XCI that is 
not critical in itself, but could reinforce the boundaries 
between active and silent chromatin along the inactive 
Xp. Because the disruption of PRC2 through the muta- 
tion of EED is associated with instability of I-XCI in 
extraembryonic derivatives [23,24], this result also im- 
plies that a function of the complex distinct from its 
H3K27me3 activity could be responsible for the main- 
tenance of Xp-linked silencing. Other PRC2 members or 
cofactors, such as JARID2, which has recently been 
shown to bridge PRC2 to Xist RNAs, may be involved 
[51,52]. PRCl-mediated repressive marks, DNA methy- 
lation or macroH2Al are other obvious candidates for 
serving this function [11,53]. Alternatively, active histone 
marks such as H3K4me2 or associated histone acetyl- 
ation may positively target XEN escapers. 

Conclusions 

Taken together, our results show that I-XCI in extraem- 
bryonic endoderm stem cells does not possess characteris- 
tics identical to X inactivation in other extraembryonic 
lineages or adult somatic tissues. Unlike the situation in 
TE cells, I-XCI in XEN cells is stably maintained in fully 
differentiated cells. We identified a large number of XEN- 
specific X inactivation escapers, which indicates that a 
mechanism dedicated to the definition of loci resistant to 
the global Xp silencing is at play in XEN cells. Import- 
antly, this mechanism does not, at least not completely, 
rely on a XEN-specific distribution of H3K27me3. The 
next goal will be to unravel the nature of this mechanism 
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and how it is recruited during initial cell commitment to 
the extraembryonic endoderm lineage. 

Methods 

Cell culture, drug treatment and ex vivo differentiation 

Male GHP7/7 and female GHP7/9 and GHP7/3 XEN cell 
lines used in this study and have been described previously 
[26]. GHP7/9 and GHP7/3 cell lines carry an Xp from a 
129.Pgkla origin and an Xm from a 129/Sv origin, provid- 
ing polymorphisms allowing for allelic distinction. XEN 
cells were grown on gelatin-coated culture dishes and rou- 
tinely cultured at 37°C in 8% C0 2 in Dulbecco's modified 
Eagle's medium (DMEM) (Life Technologies, Carlsbad, 
CA, USA) supplemented with 15% foetal bovine serum 
(Life Technologies) and 100 mM (3-mercaptoethanol [26]. 
XEN cell conversion into visceral endoderm was induced 
upon addition of 25 ng/ml BMP4 (R&D Systems, Minne- 
apolis, MN, USA) in serum containing DMEM [34]. XEN 
cell treatment with GSK126 was carried out as follows: 
three concentrations of inhibitor were tested during a 5- 
day treatment period (0.5 uM, 2 uM and 3 uM GSK126). 
Male F2 and female F3 TS cell lines and culture condi- 
tions have been described previously [12,18]. 

ChlP-chip procedure and data analysis 

ChIP assays were performed as previously described [54] 
using an H3K27me3 antibody (07-449) and an H3K4me2 
antibody (07-030; both from EMD Millipore, Billerica, 
MA, USA) and analysed as described previously [18] (see 
also Additional file 9 for the detailed protocol). ChlP-chip 
data have been deposited under Gene Expression Omni- 
bus accession number [GSE:50587]. 

Raw data were analysed using the Tobias Straub proto- 
col (http://www.epigenesys.eu/images/stories/protocols/ 
pdf/20111025114444_p43.pdf) with the Bioconductor R 
interface. This protocol includes quality assessment, data 
normalization and transformation. Raw H3K27me3 data 
in male and female livers were taken from the Gene Ex- 
pression Omnibus database [GSE:20617] [31] and sub- 
jected to the same statistical treatment as the TS and XEN 
data sets. The gene expression data of TS and XEN cells 
were extracted from [GSE:15519] [29], and gene expres- 
sion data of liver cells were extracted from [GSE:13583] 
[55]. 

RNA-FISH and Immuno-RNA-FISH 

RNA-FISH and immuno-RNA-FISH procedures were 
carried out as described at http://www.springerprotocols. 
com/Abstract/doi/10.1007/978-l-59745-406-3_18?verPrint 
=print [56] (H3K27me3 antibody, 07-449; EMD Milli- 
pore). For details of our experimental procedures, see 
Additional file 9. For probe localisation, see Additional file 
2D. Z-stacks were captured (step = 0.2 um) on a Zeiss 
Axioplan 2 microscope (Carl Zeiss Microscopy, Thornwood, 



NY, USA) equipped with a Hamamatsu ORCA ER charge- 
coupled device camera (Hamamatsu City, Japan) and con- 
trolled using Volocity Acquisition software (PerkinElmer, 
Waltham, MA, USA). Z-stacks deconvolution consisted in 
50 iterations per each channel using Huygens Professional 
software (Scientific Volume Imaging, Hilversum, the 
Netherlands). Volume measurements were performed 
using Volocity Professional software, and intensity plot 
profiles were obtained on an image analysis platform in 
ImageJ software (National Institutes of Health, Bethesda, 
MD, USA). 

Single-cell gene expression analysis 

Single-cell gene expression analysis was performed as 
described previously [18,57,58] and as recommended by 
Fluidigm (http://www.fluidigm.com/single-cell-expres- 
sion.html; South San Francisco, CA, USA). Briefly, TS 
and XEN cells were sorted by fluorescence-activated cell 
sorting using the MoFlo system (Beckman Coulter, Brea, 
CA, USA), and individual cells were distributed into 
wells or 96-well plates containing 5 ul of CellsDirect re- 
suspension buffer (Invitrogen, Carlsbad, CA, USA). The 
preamplification step consisted of 20 cycles using a mix 
of universal primer pairs to preamplify each gene simul- 
taneously. Preamplification was followed by exonuclease 
I treatment (New England Biolabs, Ipswich, MA, USA), 
and allelic qPCR was performed on a BioMark thermal 
cycler (Fluidigm). Raw efficiencies of each PCR assay 
and allelic specificity were measured on control DNA 
within each experiment. Transcript levels were extrapo- 
lated using the raw PCR efficiencies, thus allowing the 
direct comparison of different genes. Controls for allelic 
specificities of PCR assays are available upon request. 
See Additional file 2A for primer sequences. 

Additional files 



Additional file 1: Single-cell expression of lineage biomarkers in TS, 
XEN and ES female cells. (A) At 3.5 to 4.5 days postcoitum, trophoblast 
stem (TS) cells and extraembryonic endoderm (XEN) stem cells can be 
derived from the polar trophectoderm and from the primitive endoderm, 
respectively. Representative photographs of TS (F2 male cell line) and 
XEN cells (GHP7/9 female cell line) illustrating the distinct morphologies 
of the two cell types. XEN cells are highly motile and do not require 
cell-cell contact to proliferate. As expected, they exhibit two distinct 
morphologies: (1) round-shaped retractile and (2) epithelium-like [26]. 
(B) and (C) Three-dimensional projections of principal components (PC) 
(B) and heatmap (C) of single-cell steady-state RNA levels for indicated 
lineage biomarkers (16 biomarkers analysed). Hierarchical clustering shows 
significant segregation between the three cell populations (P < 1 CT 5 by 
F-test). ES cells were grown in 2i plus LIF medium [59]. Heatmap colour scale 
is the same as that shown in Figure 1 A in the main text, n = 65 female TS 
cells, n = 72 female XEN cells and n = 37 female ES cells. 

Additional file 2: (A) Primers used in single-cell gene expression 
analyses. (B) Single-cell gene expression data in undifferentiated TS and 
XEN cells. (C) Expression of, and H3K27me3 percentage enrichment and 
H3K4me2 percentage enrichment of, X-linked genes in male and female TS 
and XEN cells. (D) Genomic regions covered by the FISH probes used in the 



Merzouk ef al. Epigenetics & Chromatin 2014, 7:1 1 
http://www.epigeneticsandchromatin.eom/content/7/1/1 1 



Page 12 of 14 



study. (E) Single-cell gene expression data in TS cells differentiated for 3 or 
6 days and in XEN cells treated with BMP4 for 5 days. 

Additional file 3: H3K27me3 nuclear organisation on the inactive X 
chromosome in female TS cells. (A) Representative image of immuno- 
RNA-FISH for H3K27me3 (green) and Xist (red) on female TS cells (F3 cell 
line). Quantification of fluorescence intensities for Xist and H3K27me3 
across the inactive X domain show that the two domains do not strictly 
overlap. Maximal projections after deconvolution are shown. Scale 
bar = 5 urn. (B) Pie chart showing the percentage of nuclei exhibiting 
accumulation of Xist RNA only (red), coaccumulation of Xist RNA and 
H3K27me3 (yellow) or accumulation of H3K27me3 only (green) in female 
TS cells. (C) Boxplots showing the distribution of volumes occupied by 
Xist RNA and by H3K27me3 on the inactive X chromosome territory in 
female TS cells. The two distributions are significantly different (P< 0.05 
by Kolmogorov-Smirnov test), n > 50. 

Additional file 4: Distribution of H3K27me3 along the X-chromosomes 
in XEN and TS cells and in adult liver. (A) to (C) Scatterplots of H3K27me3 
percentages along the body of expressed X-linked genes in female 
cells/tissues (y-axis) relative to male cells/tissues (x-axis). Each dot 
represents a single gene and its respective percentage of H3K27me3 in 
the corresponding cell line or tissue, /(-means clustering was applied, 
which led to the identification of the same three classes of genes 
shown in Figure 2C in the main text: [0] = brown dots; [K27me3-high] = 
red dots and [K27me3-low] = green dots. Underneath the scatterplots, 
the pie charts show the percentage of expressed genes of each 
H3K27me3 class in each cell type. The number of expressed genes in 
each cell type is indicated (n). Liver ChlP-chip data were extracted from 
Gene Expression Omnibus ID [GSE:20617] [31] and subjected to the 
same statistical analysis as the ChlP-chip data obtained forTS and XEN 
cells. Note that only 376 X-linked genes have been analysed in the liver 
[31] compared to 642 X-linked genes in TS and XEN cell lines. We found 
a significantly higher percentage of expressed X-linked genes showing 
low levels of H3K27me3 in extraembryonic stem cells (either XEN or TS 
cells) compared to liver cells (P< 0.05 by Fisher's exact test). In contrast, 
we observed no significant difference between the two XEN cell lines 
or between either XEN cell line and TS cells (P < 0.05 by Fisher's exact 
test). Underneath the liver scatterplot, the Venn diagram shows the 
distribution of [K27-low] genes in TS cells, XEN cells and adult liver. 
This Venn diagram includes only X-linked genes that are expressed in 
all three cell types and that are common to the present study as well 
as to analyses of liver cells [31]. Expression data for TS and XEN cells 
were extracted from Gene Expression Omnibus ID [GSE:1 55 1 9] [29]. 
XEN cell lines: male GHP7/7 vs. female GHP7/3; TS cell lines: male 
F2 vs. female F3. 

Additional file 5: Distribution of H3K4me2 along the X 
chromosome in XEN cells. (A) Boxplots showing the H3K4me2 
distribution along expressed X-linked genes (Exp.) or along genes that 
are not significantly expressed (Not exp.) in male and female XEN cells 
(GHP7/7 and GHP7/9 cell lines). Similar results were obtained with female 
XEN cells (GHP7/3 cell line; not shown). Expression data were extracted 
from Gene Expression Omnibus ID [GSE:15519] [29]. n = 642 X-linked 
genes. *P<0.05 by Kolmogorov-Smirnov test. (B) Representative 
examples of H3K4me2 distribution along [K27-high] and [K27-low] genes 
in male and female XEN cells. mm9 UCSC screenshots. (C) Scatterplots of 
H3K4me2 percentages along the body of expressed X-linked genes in 
female (y-axis) relative to male (x-axis) XEN cells. Each dot represents a 
single gene and its respective percentage of H3K4me2 in the corresponding 
cell line. [K27-low] genes are shown in blue. The dotted line marks equal 
H3K4me2 percentages in male and female cells. Genes below the line are 
depleted in female compared to male cells, as expected for genes subject 
to XCI. In contrast, genes located around the line show similar levels of 
H3K4me2 in male and female cells, suggesting a biallelic enrichment in 
H3K4me2. [K27-low] genes are significantly enriched in H3K4me2 compared 
to [K27-high] genes in female XEN cells (P < 0.05 by x 2 test). 

Additional file 6: Single-cell RT-qPCR analysis on XEN cells. (A) 

Representative heatmap of a single-cell RT-qPCR BioMark chip from 
Fluidigm. Each row represents a single XEN cell (72 in total), and each 
column a specific PCR assay. This Chip shows the quantification (cycle 
threshold (Ct) values) of two different reporter genes, twelve different 



allelic assays for X-linked genes (P, paternal X-specific; M, maternal X-specific 
amplifications) and fifteen different assays for lineage-specific markers. 
Negative controls include reactions in which the reverse transcriptase 
has been omitted (RT-), reactions in the absence of cell (TE control) 
and internal controls for allele specificity of each X-linked assay performed 
on genomic DNA (gDNA) of either paternal (Pgkla) or maternal (129) 
origin. (B) Cumulative histograms showing allelic gene expression levels 
of Xist, Atp7a, Jpx, Ftx and Pbdcl on mixed populations of XEN cells 
(GHP7/9 cell line). Expression levels are standardised by the expression 
level of the reporter gene RplpO. AU, Arbitrary unit. Note that these 
levels of expression are in agreement with single-cell expression of the 
same genes shown in Figure 3 in the main text. 

Additional file 7: H3K27me3 abrogation upon treatment of female 
XEN cells with the EZH2 inhibitor GSK126. (A) Representative images 
of female XEN (GHP7/9) cells treated or not with 2 uM GSK126 for 5 days. 
No morphological differences are apparent after 5 days of treatment 
with 2 uM GSK126. (B) Western blot analysis of H3K27me3 in XEN female 
cells (GHP7/9) treated with increasing concentrations of GSK126 for 1, 3 
or 5 days. Total histone H3 detection is shown as a loading control. -, 
untreated; +, treated with GSK126. (C) RT-qPCR analysis of lineage 
markers in XEN cells (GHP7/9) treated with 2 uM GSK126 for 3 or 5 days. 
Mean + standard deviation values of the relative expression levels of 
each indicated gene are shown. Values have been standardised by the 
ubiquitously expressed RplpO gene, n = 2 independent experiments. ICM, 
Inner cell mass marker; PE, Parietal endoderm marker; PrE, Primitive 
endoderm marker; TE, Trophectoderm marker; VE, Visceral endoderm 
marker. No significant differences in the expression of most lineage 
markers were observed before and after treatment with GSK126, indicating 
that XEN cells retained their PrE identity. Only Fst was significantly 
downregulated after GSK126 treatment (P< 0.05 by x 2 test). We do 
not have any explanation for this result. (D) Representative images of 
RNA-FISH following H3K27me3 immunostaining (green) analysis in XEN 
female cells (GHP7/9 cell line) treated with GSK126 (2 uM) for 5 days. 
The inactive X chromosome is detected by Xist RNA accumulation (red). 
Primary transcription at the indicated X-linked gene is codetected in 
yellow. Scale bar = 5 um. (E) Cumulative histograms of the percentages 
of nuclei with the depicted expression pattern. Only nuclei showing a 
complete lack of H3K27me3 accumulation at the inactive X territory 
(X/sf-coated X) are scored. X-linked genes are grouped according to 
H3K27me3 level. On the left of the histogram, the diagram shows 
the localisation of RNA-FISH probes along the X chromosome. No 
significant difference in the expression pattern could be detected 
between treated and untreated cells (P< 0.05 by x 2 test), n > 50. 

Additional file 8: X-linked gene expression in differentiating female 
TS cells. (A) Representative photographs illustrating the different cell 
morphologies observed after 3 days (left panel) and after 6 days (right 
panel) of differentiation of female TS cells (F3 cell line). At day 6, TGCs (*) 
are clearly visible. (B) Heatmap of single-cell, steady-state RNA levels for 
the 20 most stringent discriminants of differentiated and undifferentiated 
TS cells (P< 10~ 3 by F-test). Hierarchical clustering has been applied. 
Complete data set is given in Additional file 2E. n = 65 female TS cells 
(green), n = 54TS differentiated for 3 days (light blue) and n = 46 female 
TS cells differentiated for 6 days (dark blue). The expected tissue specificity 
of each gene is indicated. PE, Parietal endoderm; VE, Visceral endoderm; 
ICM, Inner cell mass; EPI, Epiblast; TE, Trophectoderm; PI. Lab, Placental 
labyrinth; TGCs, Trophoblast giant cells; EPC, Fetoplacental cone. 
The colour scale is the same as that in Figure 1 A in the main text. 
(C) Three-dimensional projections of principal components (PCs) of 
single-cell expression profiles of genes shown in (B). (D) Two-colour 
RNA-FISH analysis of the indicated genes in TS cells differentiated for 
3 days. The arrowhead points to a Kdm5c transcription signal on the 
X/sf-coated inactive X chromosome that is not readily visible. Scale 
bar = 5 um. (E) Three-colour RNA-FISH analysis of the indicated genes 
in TS cells differentiated for 5 days. Arrowheads point to Rlim and 
Kdm5c transcription signals on the X/sf-coated inactive X chromosome 
that are especially difficult to see. Scale bar = 5 um. (F) Cumulative 
histogram of the percentage of nuclei with the depicted expression 
pattern in TS cells differentiated for 3 days (D3) or for 5 days (D5). 
Results obtained in undifferentiated TS cells (DO) are also shown to 
facilitate the comparison. The asterisks mark significant differences 
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between X-linked gene expression profiles in differentiated compared 
to undifferentiated TS cells (P< 0.05 by x 2 test), n > 100. 

Additional file 9: Additional Methods. 



Abbreviations 

CHART-seq: Capture hybridisation analysis of RNA targets with deep 
sequencing; ChlP-chip: Chromatin immunoprecipitation followed by chip 
hybridisation; FRAP: Fluorescence recovery after photobleaching; ICM: Inner 
cell mass; l-XCI; Imprinted X chromosome inactivation; ncRNA: Noncoding 
RNA; PE: Parietal endoderm; PrE: Primitive endoderm; RAP: RNA antisense 
purification; RNA-FISH: Fluorescence in situ hybridisation on RNA; RT-qPCR: Reverse 
transcription followed by quantitative polymerase chain reaction; 
TE: Trophectoderm; TS cell: Trophoblast stem cell; VE: Visceral endoderm; 
Xa: Active X; XEN cell: Extraendoderm stem cell; Xi: Inactive X; Xm: Maternal 
X; Xp: Paternal X. 

Competing interest 

The authors declare that they have no competing interests. 
Authors' contributions 

SM carried out some of the ChlP-chip experiments, the ChIP data analysis, 
the FISH and immunofluorescence studies, the single-cell RT-qPCR studies, 
and the XEN differentiation experiments and inhibitor treatments, and also 
participated in manuscript writing. JLD established the protocol for the 
single-cell RT-qPCR analysis, designed some of the RT-qPCR primers and 
helped to draft the manuscript. AD performed the initial cell culture and 
characterisation of XEN cell lines before the ChIP procedure. PN established 
the protocol for ChIP on XEN cells and performed the initial ChlP-chip 
experiments. PA participated in the design and coordination of the study 
and helped to write the manuscript. CM coordinated the overall experiment 
design, participated in data interpretation, designed and tested the allelic 
PCR assays, participated in immuno-RNA-FISH studies and FISH data 
analysis and wrote the manuscript. All authors read and approved the 
final manuscript. 

Acknowledgements 

We thank Michel Cohen-Tannoudji, Jerome Artus and Edith Heard for critical 
reading of the manuscript and Jerome Artus for advice on XEN cell culture 
and differentiation. This work was supported by a Marie Curie ERG grant to 
CM (224945), by recurrent funding from the Institut Pasteur and the French 
National Centre for Scientific Research (CNRS) and by grants from the French 
National Agency for Research (ANR), the Epigenome Network of Excellence, 
REVIVE Labex and the Louis D Foundation of the Institut de France (Louis D 
Prize awarded to PA). JLD was supported by a postdoctoral fellowship from 
the Association for Research against Cancer (ARC). SM was supported by 
a doctoral fellowship from the French Ministry for Scientific Research, and 
by grants from the ARC and from the REVIVE Labex. PN and AD are 
supported on a permanent basis by the Pasteur Institute. CM is supported 
on a permanent basis by the French National Institute for Scientific and 
Medical Research (INSERM). 

Author details 

'Mouse Molecular Genetics Laboratory, Pasteur Institute, 25 rue du Dr Roux, 
75015 Paris, France. 2 Pasteur Cell, Pierre and Marie Curie University (UPMC), 
25 rue du Dr Roux, 75015 Paris, France. 3 Present address: Pierre and Marie 
Curie University (UPMC), UMR7622, Institute of Biology of Paris-Seine (IBPS), 
75005 Paris, France. 4 Present address: Epigenetics of Stem Cells Laboratory', 
Pasteur Institute, 25 rue du Dr Roux, 75015 Paris, France. 5 Present address: 
Dynamics of Epigenetic Regulation, EMBL Monterotondo, Adriano 
Buzzati-Traverso Campus, Via Ramarini 32, 00015 Monterotondo, Italy. 

Received: 14 April 2014 Accepted: 2 June 2014 
Published: 20 June 2014 

References 

1. Okamoto I, Otte AP, Allis CD, Reinberg D, Heard E: Epigenetic dynamics of 
imprinted X inactivation during early mouse development. Science 2004, 

303:644-649. 



2. Okamoto I, Arnaud D, Le Baccon P, Otte AP, Disteche CM, Avner P, Heard E: 
Evidence for de novo imprinted X-chromosome inactivation independent 
of meiotic inactivation in mice. Nature 2005, 438:369-373. 

3. Patrat C, Okamoto I, Diabangouaya P, Vialon V, Le Baccon P, Chow J, Heard 
E: Dynamic changes in paternal X-chromosome activity during imprinted 
X-chromosome inactivation in mice. Proc Natl Acad Sci USA 2009, 
106:5198-5203. 

4. Kay GF, Barton SC, Surani MA, Rastan S: Imprinting and X chromosome 
counting mechanisms determine Xist expression in early mouse 
development. Cell 1 994, 77:639-650. 

5. Nesterova TB, Barton SC, Surani MA, Brockdorff N: Loss of Xist imprinting in 
diploid parthenogenetic preimplantation embryos. Dev Biol 2001, 
235:343-350. 

6. Mak W, Nesterova TB, de Napoles M, Appanah R, Yamanaka S, Otte AP, 
Brockdorff N: Reactivation of the paternal X chromosome in early mouse 
embryos. Science 2004, 303:666-669. 

7. Takagi N, Sasaki M: Preferential inactivation of the paternally derived 
X chromosome in the extraembryonic membranes of the mouse. 
Nature 1975, 256:640-642. 

8. West JD, Frels Wl, Chapman VM, Papaioannou VE: Preferential expression 
of the maternally derived X chromosome in the mouse yolk sac. 

Cell 1977, 12:873-882. 

9. Escamilla-Del-Arenal M, da Rocha ST, Heard E: Evolutionary diversity and 
developmental regulation of X-chromosome inactivation. Hum Genet 
2011, 130:307-327. 

10. Wutz A: Gene silencing in X-chromosome inactivation: advances in un- 
derstanding facultative heterochromatin formation. Nat Rev Genet 2011, 

12:542-553. 

1 1 . Brockdorff N: Noncoding RNA and Polycomb recruitment. RNA 201 3, 
19:429-442. 

12. Tanaka S, Kunath T, Hadjantonakis AK, Nagy A, Rossant J: Promotion of 
trophoblast stem cell proliferation by FGF4. Science 1998, 282:2072-2075. 

13. Calabrese JM, Sun W, Song L, Mugford JW, Williams L, Yee D, Starmer J, 
Mieczkowski P, Crawford GE, Magnuson T: Site-specific silencing of 
regulatory elements as a mechanism of X inactivation. Cell 2012, 
151:951-963. 

14. Plath K, Fang J, Mlynarczyk-Evans SK, Cao R, Worringer KA, Wang H, de la 
Cruz CC, Otte AP, Panning B, Zhang Y: Role of histone H3 lysine 27 
methylation in X inactivation. Science 2003, 300:131-135. 

15. Mak W, Baxter J, Silva J, Newall AE, Otte AP, Brockdorff N: Mitotically stable 
association of Polycomb group proteins Eed and Enxl with the inactive 
X chromosome in trophoblast stem cells. Curr Biol 2002, 12:1016-1020. 

16. Corbel C, Diabangouaya P, Gendrel AV, Chow JC, Heard E: Unusual 
chromatin status and organization of the inactive X chromosome in 
murine trophoblast giant cells. Development 2013, 140:861-872. 

1 7. Senner CE, Krueger F, Oxley D, Andrews S, Hemberger M: DNA methylation 
profiles define stem cell identity and reveal a tight embryonic- 
extraembryonic lineage boundary. Stem Cells 2012, 30:2732-2745. 

18. Dubois A, Deuve JL, Navarro P, Merzouk S, Richard S, Commere PH, Louise A, 
Arnaud D, Avner P, Morey C: Spontaneous reactivation of clusters of 
X-linked genes is associated with the plasticity of X-inactivation in 
mouse trophoblast stem cells. Stem Cells 2014, 32:377-390. 

19. Hadjantonakis AK, Cox LL, Tarn PP, Nagy A: An X-linked GFP transgene 
reveals unexpected paternal X-chromosome activity in trophoblastic 
giant cells of the mouse placenta. Genesis 2001, 29:133-140. 

20. Hemberger M, Kurz H, Orth A, Otto S, Luttges A, Elliott R, Nagy A, Tan SS, 
Tarn P, Zechner U, Fundele RH: Genetic and developmental analysis of 
X-inactivation in interspecific hybrid mice suggests a role for the Y 
chromosome in placental dysplasia. Genetics 2001, 157:341-348. 

21 . Tarn PP, Williams EA, Tan SS: Expression of an X-linked HMG-lacZ 
transgene in mouse embryos: implication of chromosomal imprinting 
and lineage-specific X-chromosome activity. Dev Genet 1994, 15:491-503. 

22. Tan SS, Williams EA, Tarn PP: X-chromosome inactivation occurs at different 
times in different tissues of the post-implantation mouse embryo. 
Nat Genet 1993, 3:170-174. 

23. Kalantry S, Mills KC, Yee D, Otte AP, Panning B, Magnuson T: The Polycomb 
group protein Eed protects the inactive X-chromosome from 
differentiation-induced reactivation. Nat Cell Biol 2006, 8:195-202. 

24. Wang J, Mager J, Chen Y, Schneider E, Cross JC, Nagy A, Magnuson T: 
Imprinted X inactivation maintained by a mouse Polycomb group gene. 
Nat Genet 2001,28:371-375. 



Merzouk ef al. Epigenetics & Chromatin 2014, 7:1 1 
http://www.epigeneticsandchromatin.eom/content/7/1/1 1 



Page 14 of 14 



25. Fowler KJ, Mitrangas K, Dziadek M: In vitro production of Reichert's 
membrane by mouse embryo-derived parietal endoderm cell lines. 
Exp Cell Res 1990, 191:194-203. 

26. Kunath T, Arnaud D, Uy GD, Okamoto I, Chureau C, Yamanaka Y, Heard E, 
Gardner RL, Avner P, Rossant J: Imprinted X-inactivation in extra- 
embryonic endoderm cell lines from mouse blastocysts. Development 
2005, 132:1649-1661. 

27. Niakan KK, Schrode N, Cho LT, Hadjantonakis AK: Derivation of 
extraembryonic endoderm stem (XEN) cells from mouse embryos and 
embryonic stem cells. Nat Protoc 2013, 8:1028-1041. 

28. Ohhata T, Senner CE, Hemberger M, Wutz A: Lineage-specific function of 
the noncoding Tsix RNA for Xist repression and Xi reactivation in mice. 
Genes Dev 2011, 25:1702-1715. 

29. Rugg-Gunn PJ, Cox BJ, Ralston A, Rossant J: Distinct histone modifications 
in stem cell lines and tissue lineages from the early mouse embryo. 
Proc Natl Acad Sci USA 2010, 1 07:1 0783-1 0790. 

30. Greenfield A, Carrel L, Pennisi D, Philippe C, Quaderi N, Siggers P, Steiner K, 
Tarn PP, Monaco AP, Willard HF, Koopman P: The UTXgene escapes X 
inactivation in mice and humans. Hum Mol Genet 1998, 7:737-742. 

31. Yang F, BabakT, Shendure J, Disteche CM: Global survey of escape from X 
inactivation by RNA-sequencing in mouse. Genome Res 2010, 20:614-622. 

32. McCabe MT, Ott HM, Ganji G, Korenchuk S, Thompson C, Van Aller GS, Liu Y, 
Graves AP, Delia Pietra A 3rd, Diaz E, LaFrance LV, Mellinger M, Duquenne C, 
Tian X, Kruger RG, McHugh CF, Brandt M, Miller WH, Dhanak D, Verma SK, 
Tummino PJ, Creasy CL: EZH2 inhibition as a therapeutic strategy for 
lymphoma with EZH2-activating mutations. Nature 2012, 492:108-1 12. 

33. Sado T, Fenner MH, Tan SS, Tam P, Shioda T, Li E: X inactivation in the 
mouse embryo deficient for Dnmtl: distinct effect of hypomethylation 
on imprinted and random X inactivation. Dev Biol 2000, 225:294-303. 

34. Artus J, Douvaras P, Piliszek A, Isern J, Baron MH, Hadjantonakis AK: 
BMP4 signaling directs primitive endoderm-derived XEN cells to an 
extraembryonic visceral endoderm identity. Dev Biol 2012, 361:245-262. 

35. Paca A, Seguin CA, Clements M, Ryczko M, Rossant J, Rodriguez TA, Kunath 
T: BMP signaling induces visceral endoderm differentiation of XEN cells 
and parietal endoderm. Dev Biol 2012, 361:90-102. 

36. Yang C, McLeod AJ, Cotton AM, de Leeuw CN, Laprise S, Banks KG, 
Simpson EM, Brown CJ: Targeting of >1.5 Mb of human DNA into the 
mouse X chromosome reveals presence of cis-acting regulators of 
epigenetic silencing. Genetics 2012, 192:1281-1293. 

37. Li N, Carrel L: Escape from X chromosome inactivation is an intrinsic 
property of the Jarid 1c locus. Proc Natl Acad Sci USA 2008, 1 05:1 7055-1 7060. 

38. Horvath LM, Li N, Carrel L: Deletion of an X-inactivation boundary disrupts 
adjacent gene silencing. PLoS Genet 2013, 9:e1 003952. 

39. Filippova GN, Cheng MK, Moore JM, Truong JP, Hu YJ, Nguyen DK, Tsuchiya 
KD, Disteche CM: Boundaries between chromosomal domains of X 
inactivation and escape bind CTCF and lack CpG methylation during 
early development. Dev Cell 2005, 8:31-42. 

40. Lopes AM, Arnold-Croop SE, Amorim A, Carrel L: Clustered transcripts that 
escape X inactivation at mouse XqD. Mamm Genome 201 1, 22:572-582. 

41. Tsuchiya KD, Greally JM, Yi Y, Noel KP, Truong JP, Disteche CM: 
Comparative sequence and X-inactivation analyses of a domain of 
escape in human Xpll.2 and the conserved segment in mouse. 
Genome Res 2004, 14:1275-1284. 

42. Carrel L, Park C, Tyekucheva S, Dunn J, Chiaromonte F, Makova KD: 
Genomic environment predicts expression patterns on the human 
inactive X chromosome. PLoS Genet 2006, 2:e15l. 

43. Reinius B, Shi C, Hengshuo L, Sandhu KS, Radomska KJ, Rosen GD, Lu L, 
Kullander K, Williams RW, Jazin E: Female-biased expression of long 
non-coding RNAs in domains that escape X-inactivation in mouse. 
BMC Genomics 2010, 11:614. 

44. Ciavatta D, Kalantry S, Magnuson T, Smithies O: A DNA insulator prevents 
repression of a targeted X-linked transgene but not its random or 
imprinted X inactivation. Proc Natl Acad Sci USA 2006, 103:9958-9963. 

45. Berletch JB, Yang F, Xu J, Carrel L, Disteche CM: Genes that escape from 
X inactivation. Hum Genet 201 1, 130:237-245. 

46. Carrel L, Willard HF: X-inactivation profile reveals extensive variability in 
X-linked gene expression in females. Nature 2005, 434:400-404. 

47. Lingenfelter PA, Adler DA, Poslinski D, Thomas S, Elliott RW, Chapman VM, 
Disteche CM: Escape from X inactivation of Smcx is preceded by 
silencing during mouse development. Nat Genet 1998, 18:212-213. 



48. Prothero KE, Stahl JM, Carrel L: Dosage compensation and gene 
expression on the mammalian X chromosome: one plus one does not 
always equal two. Chromosome Res 2009, 17:637-648. 

49. Marks H, Chow JC, Denissov S, Francoijs KJ, Brockdorff N, Heard E, 
Stunnenberg HG: High-resolution analysis of epigenetic changes 
associated with X inactivation. Genome Res 2009, 19:1361-1373. 

50. Cerase A, Smeets D, Tang YA, Gdula M, Kraus F, Spivakov M, Moindrot B, 
Leleu M, Tattermusch A, Demmerle J, Nesterova TB, Green C, Otte AP, 
Schermelleh L, Brockdorff N: Spatial separation of Xist RNA and polycomb 
proteins revealed by superresolution microscopy. Proc Natl Acad Sci USA 
2014, 111:2235-2240. 

51 . Casanova M, Preissner T, Cerase A, Root R, Yamada D, Li X, Appanah R, 
Bezstarosti K, Demmers J, Koseki H, Brockdorff N: Polycomblike 2 facilitates 
the recruitment of PRC2 Polycomb group complexes to the inactive X 
chromosome and to target loci in embryonic stem cells. Development 

2011, 138:1471-1482. 

52. da Rocha ST, Boeva V, Escamilla-Del-Arenal M, Ancelin K, Granier C, Matias NR, 
Sanulli S, Chow J, Schulz E, Picard C, Kaneko S, Helin K, Reinberg D, Stewart AF, 
Wutz A, Margueron R, Heard E: jarid2 is implicated in the initial Xist-induced 
targeting of PRC2 to the inactive X chromosome. Mol Cell 2014, 53:301-316. 

53. Changolkar LN, Singh G, Cui K, Berletch JB, Zhao K, Disteche CM, Pehrson 
JR: Genome-wide distribution of macroH2A1 histone variants in mouse 
liver chromatin. Mol Cell Biol 2010, 30:5473-5483. 

54. Navarro P, Richard S, Ciaudo C, Avner P, Rougeulle C: Tsix transcription 
across the Xist gene alters chromatin conformation without affecting Xist 
transcription: implications for X-chromosome inactivation. Genes Dev 
2005, 19:1474-1484. 

55. Fujimoto T, Miyasaka K, Koyanagi M, Tsunoda T, Baba I, Doi K, Ohta M, Kato 
N, Sasazuki T, Shirasawa S: Altered energy homeostasis and resistance to 
diet-induced obesity in Kft/tP-deficient mice. PLoS One 2009, 4:e4240. 

56. Chaumeil J, Okamoto I, Heard E: X-chromosome inactivation in mouse 
embryonic stem cells: analysis of histone modifications and 
transcriptional activity using immunofluorescence and FISH. Methods 
Enzymol 2004, 376:405-419. 

57. Guo G, Huss M, Tong GQ, Wang C, Sun LL, Clarke ND, Robson P: Resolution 
of cell fate decisions revealed by single-cell gene expression analysis 
from zygote to blastocyst. Dev Cell 2010, 18:675-685. 

58. Rugg-Gunn PJ, Cox BJ, Lanner F, Sharma P, Ignatchenko V, McDonald AC, 
Garner J, Gramolini AO, Rossant J, Kislinger T: Cell-surface proteomics 
identifies lineage-specific markers of embryo-derived stem cells. Dev Cell 

2012, 22:887-901. 

59. Silva J, Barrandon O, Nichols J, Kawaguchi J, Theunissen TW, Smith A: 
Promotion of reprogramming to ground state pluripotency by signal 
inhibition. PLoS Biol 2008, 6:e253. 



doi:1 0.1 186/1 756-8935-7-11 

Cite this article as: Merzouk ef al:. Lineage-specific regulation of 
imprinted X inactivation in extraembryonic endoderm stem cells. 

Epigenetics & Chromatin 2014 7:1 1 . 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www. biomedcentra I .com/su bmit 



o 



BioMed Central 



